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LONG-TERM  GOALS 

(a)  To  provide  an  aeeurate,  effieient,  eomputational  model  (two-seale  approximation,  TSA)  for  the 
4-wave  interactions,  in  operational  wave  forecast  models,  suitable  for  global,  basin  and  coastal 
scale  applications,  and  able  to  transition  seamlessly  from  deep  to  shallow  water. 

(b)  Fully  test  TSA  with  respect  to  exact  codes  for  the  full  Boltzmann  integral  (FBI),  for  duration- 
limited,  fetch-limited  wave  growth,  turning  winds,  swell-windsea,  interactions,  etc. 

(c)  Numerically  investigate  and  clarify  the  basis  for  TSA,  its  limitations,  errors,  enhancements, 
improvements,  self-similarity  properties,  and  spectral  flux  properties. 

(d)  Implement  TSA  in  a  variety  of  modem  operational  wave  forecast  models,  e.g.  WAVEWATCH™ 
(WW3)  and  SWAN  for  extensive  tests  on  important,  realistic  wave  conditions. 

(e)  Derive,  adapt  and  implement  new  formulations  for  source  terms.  Sin  and  Sds,  from  recent 
literature  and  the  NOPP  partnership,  with  TSA,  in  modern  wave  models,  for  tests,  including 
veering  or  accelerating  winds,  sea  and  swell  interactions,  and  real  storm  cases. 

OBJECTIVES 

For  this  reporting  period: 

1)  To  formulate  a  relatively  extensive  set  of  reference  spectra  with  corresponding  FBI  results,  in  order 
to  facilitate  definitive  comparisons  with  TSA,  and  to  lead  to  better  understanding  of  limitations  and 
areas  in  need  of  improvement  in  the  TSA  formulation. 

2)  Results  for  TSA  (in  comparison  with  FBI)  for  directionally  varying  spectra,  rapidly  veering  wind 
conditions  and  other  difficult  spectral  cases  where  the  formulation  may  fail. 

3)  Preliminary  implementation  of  TSA  in  WW3,  re-stmcturing  of  code  as  necessary 

APPROACH 

The  project  initially  needs  to  focus  on  TSA’s  performance  with  respect  to  the  fundamental 
properties  of  the  FBI  (full  Boltzmann  integral)  code. 

1)  Reference  spectra  and  FBI,  The  spectral  tests  presented  in  Resio  and  Perrie  (2008)  and  Perrie  and 
Resio  (2009)  are  limited  and  somewhat  idealized.  Further  work  is  needed  to  develop  a  reasonable 
set  of  reference  spectra  and  their  corresponding  FBI  results,  in  order  to  facilitate  more  definitive 
comparisons  for  TSA.  This  would  lead  the  way,  in  terms  of  clarifying  numerical  limitations  of 
the  present  TSA  formulation,  and  determination  of  areas  where  further  development  is  needed. 
This  activity  is  done  by  Resio  (USAGE),  Eong  (USAGE)  and  Perrie  (BIO). 
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2)  Reference  spectra  for  TSA,  How  does  TSA  perform  for  direetionally  varying  (sheared)  speetra, 

for  example  in  rapidly  veering  winds,  eompared  to  fundamental  properties  of  the  FBI  eode? 

Perrie  and  Resio  (2009)  show  that  this  test  is  eritieal  for  TSA,  using  buoy  data  from  hurrieane 
Wilma.  Clearly  TSA  has  ability  to  reproduee  results  suggested  by  FBI,  but  in  more  extreme  eases 
shows  some  evidenee  of  divergenee  from  FBI  in  eritieal  tests.  We  will  test  diffieult  speetral 
eases,  in  order  to  elarify  limitations  of  the  present  TSA  formulation,  and  to  determine  areas 
where  more  TSA  development  and  improvements  are  needed.  This  work  is  done  by  Resio 
(USAGE),  Long  (USAGE)  and  Perrie  (BIO). 

3)  TSA  put  into  WW3.  We  will  make  a  preliminary  implementation  of  TSA  in  the  latest  WW3 

version.  This  is  relatively  straightforward  beeause  we  have  used  versions  of  WW3  for  many 
years  in  the  Gulf  of  Maine  Oeean  Observing  System  -  GoMOOS;  moreover,  it  is  preliminary  in 
the  sense  that  TSA  is  going  to  evolve  and  be  enhaneed,  with  bugfixes  implemented,  as  required, 
beeoming  more  optimal  as  the  projeet  progresses.  This  task  is  preparatory  in  laying  the  basis 
needed  to  implement  TSA  in  operational  wave  models;  and  yet  this  work  is  fundamental, 
requiring  re-strueturing  of  the  TSA  eode  to  eomply  with  WW3  arehiteeture.  This  work  is  done  by 
Toulany  (BIO),  Auelair  (BIO/  Dalhousie),  Resio  (USAGE),  Long  (USAGE),  Perrie  (BIO),. 

WORK  COMPLETED 

Eunding  delays  have  oeeurred  beeause  polieies  were  unelear  regarding  eollaborative  agreements  that 

would  be  needed  between  Bedford  Institute  of  Oeeanography  (Government  of  Canada)  and  ONR. 

These  diffieulties  are  now  largely  resolved.  The  projeet  will  move  ahead  on  sehedule,  beeause  work 

was  eompleted  using  leveraged  funds  from  a  eomplementary  Canadian  Government  projeet  funded  by 

the  Panel  on  Energy  Researeh  and  Development  whieh  is  ending  on  3 1  Mareh  2011. 

1)  Reference  spectra  and  FBI,  The  pioneering  work  of  Zakharov  and  Eilonenko  (1966)  showed  that 
nonlinear  wave-wave  interaetions  produee  stationary  solutions  related  to  eonstant  energy  flux 
through  the  equilibrium  range  when  a  deep-water  speetrum  follows  an  form,  as  verified  in 
numerieal  studies  in  whieh  the  direetional  speetra  follows  a  eonstant  angular  spreading  distribution. 
In  Resio  et  al.  (2010)  we  show  that,  although  energy  fluxes  through  sueh  speetra  remain  essentially 
eonstant,  momentum  fluxes  do  not.  Moreover,  if  the  angular  distribution  of  the  speetrum  is  allowed 
to  behave  in  a  manner  eonsistent  with  observations,  both  the  energy  fluxes  and  the  momentum 
fluxes  appear  to  remain  eonstant  through  a  major  portion  of  the  speetrum.  Thus,  it  appears  that  the 
direetional  distribution  of  energy  within  wind  wave  speetra  has  adjusted  to  a  form  with  non- 
divergent  nonlinear  fluxes,  suggesting  that  these  fluxes  likely  play  a  very  prominent  role  in  the 
evolution  of  direetional  speetra  during  wave  generation. 

2)  Reference  spectra  for  TSA,  A  preliminary  set  of  eomputations  were  eompleted  showing  how 
TSA  performs  for  direetionally  varying  (sheared)  speetra,  in  rapidly  veering  winds,  and  in  swell- 
wave  interaetion  situations,  eompared  to  fundamental  properties  of  the  FBI  eode.  This  aetivity 
extends  basie  eases  eonsidered  by  Perrie  and  Resio  (2009)  in  using  observed  wave  data  from 
hurrieane  Wilma.  We  are  able  to  show  that  TSA  has  ability  to  reproduee  basie  results  suggested  by 
FBI,  even  in  very  extreme  eases.  We  will  eontinue  to  test  different  speetral  eases,  in  order  to 
determine  areas  where  TSA  improvements  are  needed. 

3)  TSA  put  into  WW3,  We  have  made  a  preliminary  implementation  of  TSA  in  the  latest  WW3 
version  (3.14).  This  aetivity  has  involved  eonstruetion  of  the  diagonal  term  (WAMDIG,  1988) 
needed  for  the  implieit  integration  used  by  WW3.  Preliminary  results  are  shown  in  Figure  1  below. 
These  results  are  eneouraging,  in  that  we  ean  generate  wave  grow  eurves  that  aet  reasonably. 
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although  the  other  souree  terms  for  wind  input  and  dissipation  have  not  been  tuned  to  TSA. 
Additional  work  needs  to  address  outstanding  issues. 

RESULTS 

a)  In  Resio  et  al.  (2010)  we  examine  nonlinear  fluxes  of  energy  and  momentum  through  wave  speetra 
via  an  exaet  integration  of  the  full  Boltzmann  integral  (the  WRT  method).  Conelusions  are: 

i.  The  bimodal  strueture  observed  in  studies  of  the  direetional  distributions  of  wind-wave  speetra 
(Wang  and  Hwang,  2001;  Long  and  Resio,  2007;  Toffoli  et  al,  2010)  is  relatively  eonsistent  with 
the  eos^"  angular  distributions,  at  least  in  a  bulk  sense,  as  derived  in  earlier  field  studies. 

ii.  Speetra  with  a  bimodal  distribution  of  energy  eonsistent  with  reeent  observations  give  relatively 
eonstant  fluxes  of  both  energy  and  momentum  through  the  equilibrium  range  whieh  suggests  that 
the  role  of  nonlinear  interaetions  is  quite  eritieal  to  the  direetional  evolution  of  wave  speetral. 

ill.  Nonlinear  momentum  fluxes  from  the  speetral  peak  region  through  the  equilibrium  range,  from  a 
numerieal  solution  to  the  Boltzmann  integral,  agree  well  with  the  expeeted  momentum  balanee. 

b)  In  studies  of  referenee  speetra  and  TSA,  we  showed  that  for  direetionally  varying  (sheared)  speetra, 
in  rapidly  veering  winds,  and  in  swell-wave  interaetion  situations,  eompared  to  fundamental 
properties  of  the  FBI  eode,  TSA  has  ability  to  reproduee  basie  results  suggested  by  FBI,  even  in 
very  extreme  eases.  We  will  eontinue  to  test  different  speetral  eases. 

IMPACT/APPLICATIONS 

This  projeet  is  eoneerned  with  TSA.  The  rationale  for  this  projeet  is  that  nonlinear  4-wave  interaetions, 
represented  by  FBI,  have  been  shown  to  be  eentral  to  wave  foreeast  models  sinee  Hasselmann  (1962) 
and  Zakharov  and  Filonenko  (1966),  and  aeeurate  numerieal  formulations  have  been  demonstrated,  by 
Resio  and  Perrie  (2008),  and  Resio  et  al  (2001),  for  finite  depths.  However,  these  formulations  for  the 
FBI  are  too  slow  for  operational  wave  foreeasting.  The  diserete  interaetion  approximation,  or  DIA 
(WAMDIG,  1988)  is  the  basis  for  operational  wave  models,  sueh  as  WW3  and  SWAN.  Although  DIA 
is  ealibrated  to  agree  with  the  FBI  for  the  forward  faee  of  the  speetrum,  e.g.  with  peakedness  y  =  3.3,  it 
deviates  substantially  elsewhere  within  the  speetrum.  On  the  other  hand  TSA  offers  a  new  approaeh  to 
this  problem.  TSA  is  presented  by  Resio  and  Perrie  (2008),  and  Perrie  and  Resio  (2009).  It 
outperforms  DIA  in  aeouraey,  with  respeet  to  exaet  FBI  results,  and  it  has  run  times  that  are  on  the 
seale  of  run  times  for  DIA;  therefore  TSA  is  a  eandidate  for  implementation  in  WW3  and  SWAN. 

RELATED  PROJECTS 

A  related  projeet  is  funded  by  the  Canadian  Panel  on  Energy  Researeh  and  Development  entitled 
“Waves  and  Winds  in  Extreme  Storms”.  It’s  foeus  is  development  of  a)  new  wave  model  physies,  b) 
improved  wind  and  waves  estimates  from  better  atmosphere-oeean  eoupling,  e)  improved  models  for 
wave-eurrent  interaetions,  d)  development  of  a  versatile  prototype  wave  foreeast  system,  e)  estimates 
of  biases  in  wind  and  wave  elimatologies.  This  projeet  is  ending  on  31  Mareh  2011. 
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Eigure  1.  Variation  in  (a)  significant  wave  height  Hs,  and  (b)  peak  frequency  fp,  as  a  function  of  time  t, 
for  an  integration  of  TSA  implemented  in  WW3  version  3.14,  implemented  in  a  square-ocean 
SWAMP  test,  compared  to  results  using  DIA  and  WRT  for  the  nonlinear  interactions. 


4 


